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A power-law description of heavy ion collision centrality 



O 

o 

C 

a 

o 

m 

m 
> 



O , 
^ : 

i 

^ ■ 



X 



Thomas A. Trainor and Duncan J. Prindlc 

CENPA, Box 354290, University of Washington, Seattle, WA 98195 
(Dated: February 2, 2008) 

The minimum-bias distribution on heavy ion collision multiplicity n cn is well approximated by 
power- law form n~^*, suggesting that a change of variable to nj£ may provide more precise access 
to the structure of the distribution and to A-A collision centrality. We present a detailed centrality 
study of Hijing-1.37 Monte Carlo data at 200 GeV using the power-law format. We find that the 



minimum-bias distribution on n 



1/4 



determined with a Glauber Monte Carlo simulation, is 



participant ' 

uniform except for a 5% sinusoidal variation. The power-law format reveals precise linear relations 
between Glauber parameters n par t and nun and the fractional cross section. The power-law format 
applied to RHIC data facilitates incorporation of extrapolation constraints on data and Glauber 
distributions to obtain a ten-fold improvement in centrality accuracy for peripheral collisions. 

PACS numbers: 24.60.Ky,25.75.Gz 

Keywords: Collision centrality, heavy ion collisions, participant scaling, binary-collision scaling, Glauber 
model, two-component model, Hijing simulation 



I. INTRODUCTION 



Two-particle correlations in RHIC heavy ion collisions 
change rapidly with collision centrality, reflecting strong 
changes in collision dynamics [jj, [2|, y, |J] . Some variations 
relate to the changing geometry of hadronization [3J. 
Other variations arise from copious low-Q 2 parton scat- 
tering [5|. Correlations from p-p collisions with simpler 
dynamics provide a precision reference for peripheral Au- 
Au collisions Q. However, conventional RHIC centrality 
methods have been ineffective for the 20% most periph- 
eral collisions, where changes relative to p-p are rapid 
and informative. We are motivated therefore to improve 
centrality determination for A-A collisions. 



The power-law representation of the minimum-bias dis- 
tribution on collision multiplicity provides a basis for ma- 
jor improvement. The minimum-bias distribution, con- 
ventionally plotted in a semi- log format, is approximately 
a power-law distribution 7], implying that a change of 
variable should lead to a more compact form of the dis- 
tribution. The ideal power-law form is a uniform distri- 
bution between two well-defined endpoints. Deviations 
from that precision reference are easily identified and 
studied. 



The power-law format helps in several ways. 1) It pro- 
vides accurate centrality determination down to N-N col- 
lisions, even if the measured 'minimum-bias' distribution 
is strongly distorted or biased by triggering and vertex 
reconstruction inefficiencies. 2) Applied to the Glauber 
model the power-law format provides compact and pre- 
cise representations of the Glauber parameters. 3) The 
power-law plotting format reveals distribution details at 
the few-percent level important for precision comparisons 
with p-p collisions. 



II. METHOD 

The novel techniques described in this paper emerged 
from the observation that minimum-bias distributions on 
n P art (participant nucleon number), n c h, Pt and E t vary 
approximately as x -3 / 4 , a power law. Transforming the 
distributions to x 1 ^ 4 has led to major improvements in 
analysis accuracy. Through running integrals of several 
power-law distributions we are able to relate measured 
quantities and geometry parameters at the percent level. 
The methods are applied to Hijing data and RHIC data. 

We first introduce the power-law form of the minimum- 
bias distribution, compare it to the conventional semi-log 
form and describe its properties. We then review the 
Glauber model of A-A collision geometry and study the 
power-law form of the Glauber minimum-bias distribu- 
tions on ripart and n^ n . We construct parameterizations 
of running integrals which relate the Glauber parame- 
ters to the fractional cross section. We then demon- 
strate power-law centrality determination with Hijing 
and RHIC data, relating n c h to fractional cross section 
with percent errors. 

To demonstrate the overall method we study the cen- 
trality dependence of particle and pt production in Hi- 
jing for two event classes (quench-on, quench-off) and 
particle and E t production in RHIC data. Such produc- 
tion studies provide the most demanding test of centrality 
precision by comparing centralities from data with cen- 
tralities from a Glauber Monte Carlo. Finally, we make 
a detailed comparison of systematic centrality errors in 
conventional and power-law contexts. 

We also include three appendices: A) Particle- 
production algebra: a unified approach to particle pro- 
duction in different plotting contexts; B) Numerical inte- 
gration techniques: methods of binning, numerical inte- 
gration and running integrals; C) Relating centrality pa- 
rameters: what lies behind centrality methods - joint and 
marginal distributions, running integrals of marginals, 
the participant-scaling model and the role of fluctuations. 
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III. HIJING DATA 

Hijing-1.37 § was used to produce minimum-bias 
event ensembles with ~ 1M total events for each of two 
classes: 1) quench-off Hijing - jet production but no jet 
quenching and 2) quench-on Hijing jet production with 
jet quenching. Charged particles with pseudorapidity 
\rj\ < 1, transverse momentum p t E [0.15,2] GeV/c and 
full azimuth were accepted. The total charged-particle 
multiplicity and transverse momentum in the acceptance 
for each event defines centrality parameters n c h and pt- 
The variations with centrality of Hijing quench-on and 
quench-off event classes are distinguishable, but to an 
extent depending on the analysis method. We use the 
two Hijing classes to illustrate the sensitivity of different 
centrality formats to particle and pt production mecha- 
nisms, as well as the general features of power-law cen- 
trality determination. 



IV. THE CONVENTIONAL MINIMUM-BIAS 
DISTRIBUTION 

In Fig. [1] (left panel) we plot the minimum-bias dis- 
tribution of event number on multiplicity n c h in a con- 
ventional semi-log format for the two Hijing event types. 
The distributions are monotonically decreasing, and the 
tail widths at large n c h reflect fluctuations in particle 
production for central collisions (6 = 0). The fluctuation 
magnitude (relative to Poisson) depends on the detector 
acceptance. Comparing the two Hijing event types, the 
only apparent difference is the end-point positions. 




FIG. 1: Left panel: Minimum-bias distributions for quench- 
off and quench-on (larger endpoint) Hijing collisions plotted 
with a conventional semi-log plotting format. Right panel: 
The same distributions plotted in a log- log format. The line 
represents an n~^ A power-law trend. 

In Fig. Q] (right panel) the same data are plotted in 
a log-log format which reveals that the data are well 
approximated by power-law trend n~^ 4 (solid line) 0. 
The power-law trend suggests that a change of variable 
could provide more precise access to distribution struc- 
ture. If da/dn c h oc n c ^ 4 then n^ h 4 da / dn c h ~ constant. 
Interpreting n^ 4 as the Jacobian of a variable transfor- 



mation, i.e., dn c h = 4:nJ h dnj h , we expect da/dnj h ~ 
constant. We convert the minimum-bias Hijing data dis- 
tributions to power-law form da/dnj^ = 4n^ 4 da / dn c h 
and thereby obtain more precise access to data and im- 
proved centrality determination. 



V. THE POWER-LAW DISTRIBUTION 

We introduce the power-law minimum-bias distribu- 
tion and consider examples from Hijing Monte Carlo data 
and RHIC data. We identify the principal features of the 
distribution in comparison to a participant-scaling refer- 
ence. We parameterize the form of the distribution based 
on a two-component model of particle production. 



A. Hijing 

Fig- El (left panel) shows distributions da/dnj^ vs nj^ 
for two Hijing configurations which confirm the basic fea- 
tures anticipated for the power-law format — an approx- 
imately rectangular distribution with limited amplitude 
variation and well-defined endpoints. The distribution on 
integers has been rebinned to 50 uniform bins on nj } 4 (cf. 
App. |B~3|) . insuring nearly-uniform statistical bin errors 
while retaining adequate resolution near the upper end- 
point. The solid points at the lower endpoint indicate the 
limiting edge resolution. The first few points are defined 
by the smallest integers, not the rebinning. The lower 
endpoint (half- maximum point) is n p ~ unn/2 (but cf. 
App. IC3[) . The upper endpoints n =1210 and 1400 
are not exactly at the half-maximum points due to the 
asymmetric (skewed) shape of fluctuations for central col- 
lisions as modeled by Hijing. The sloped dash-dot lines 
starting at and approximating the data are defined 
in the next subsection. 

The five parameters which describe the shape of the 
power-law minimum-bias distribution are summarized in 
Fig. [2] (right panel): 1) Upper half-maximum endpoint 
no estimates the mean n c h corresponding to central A-A 
collisions (b = 0) and upper endpoint n part .a/2 of the 
minimum-bias distribution on participant-pair number 
n par t/2 (cf. Sec. IVI CI for n par t endpoint definitions). 
2) Lower half-maximum endpoint n p is approximately 
one-half the mean n c h for N-N (~ p-p) collisions (lower 
endpoint n part , p /2 ~ 1/2). 3), 4) The slopes at the end- 
points measure particle-production fluctuations in N-N 
and 6 = A-A collisions, which also depend on the de- 
tector acceptance [13]. 5) The slope near the midpoint 
reflects the particle-production mechanism, i.e., the rela- 
tive importance of binary-collision and participant-pair 
scaling. At 200 GeV the Hijing (Pythia) N-N multi- 
plicity is njvAf ~ 5 in \rj\ < 1. For a symmetric N-N 
multiplicity distribution the lower half-maximum point 
would be n p = njviv/2. However, for Hijing (and data) 
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FIG. 2: Left panel: The power-law form of the minimum- 
bias distribution for quench-on (solid curve) and quench-off 
(dash-dot curve) Hijing data. The solid rectangles are power- 
law references for the data. The dashed rectangle represents 
participant scaling. The dash-dot lines are as defined for 
Fig. |3] (right panel) . Right panel: A cartoon of the power-law 
minimum-bias distribution identifying the principal features. 
The vertical solid lines indicate the positions of half-maximum 
endpoints n^ 4 and n^ 4 . The dashed line indicates the posi- 
tion for the N-N multiplicity njyjv ~ 2n p . 



was corrected for trigger and tracking inefficiencies and 
backgrounds. Its analysis illustrates some aspects of the 
collision geometry problem relating to real data. The 
left panel shows the conventional semi-log plotting for- 
mat with uniform bins on n h - . The assumed total cross 
section is ctq — 7.2 barns Mi. The event-trigger efficiency 
(coincidence of two ZDCs [12|) was greater than 98% for 
all multiplicities. 

The measured event-vertex efficiency was 100% for 
n h - > 50 but dropped to 60% below n h - = 5. The 
contribution to do from the lowest bin was estimated to 
be 21% based on an extrapolation using Hijing. It was 
argued that peripheral A-A collisions are linear super- 
positions of N-N collisions, and Hijing (Pythia) models 
N-N collisions correctly. Hijing was normalized to data 
for n h - G [5,25] (ji^ 4 S [1.5,2.24]) and used for the ex- 
trapolation below n h - = 5. The estimated systematic 
uncertainty in the inferred differential cross section was 
10%, due to uncertainties in ctq and the inferred rela- 
tive contribution from the first bin based on the Hijing 
extrapolation. 



n p ~ 0.45 n^w because the N-N distribution is signifi- 
cantly skewed (0.45 estimates n partiP /2, cf. App. IC 3)) . 

The power-law reference distribution is a rectangle 
(uniform distribution on nY 4 ) bounded by endpoints 
Up 4 and n\j 4 , with area ctq the total cross section defined 



The average value of the power-law 

V4_^l/4^ 
P 



by the event trigger 

minimum-bias distribution is therefore ao/{n' , ' i — rip*) 
The solid rectangles in Fig. [5] (left panel) are power- 
law references for Hijing data, with upper endpoints 
uq = 1210 (1400) for quench-off (quench-on) events 
within two units of pseudorapidity. The dashed rectan- 
gle is a participant scaling reference with lower endpoint 
n p = n par t :P /2 ■ riNN (in common with the solid rectan- 



gles) and upper endpoint no 



lpart,o/2 • flNN- 



Systematic deviations from the ideal power-law distri- 
bution near n p could result from trigger and/or vertex 
inefficiencies, or contamination from non-hadronic back- 
grounds. The five shape parameters provide precise de- 
termination of collision centrality and particle (or p t , E t ) 
production. More detailed features are suppressed in the 
running integrals used to relate n c h, Pt and E t to frac- 
tional cross section er/ero- 

Fig. [5] (left panel) demonstrates that the power-law 
format clearly reveals physics-related differences be- 
tween quench-on and quench-off distributions at the few- 
percent level not apparent in the semi-log plotting for- 
mat of Fig. [T] We now apply the power-law format to a 
minimum-bias distribution obtained from RHIC data. 



B. RHIC data 

Fig.[3]shows a minimum-bias distribution on n h - (neg- 
ative hadron multiplicity) in \rj\ < 0.5 for 60k Au-Au 
collisions at ^/s^n — 130 GeV :12j. The distribution 




FIG. 3: Left panel: A minimum-bias distribution of RHIC 
data with a conventional semi-log plotting format. Right 
panel: The same data plotted with a power-law format. The 
solid rectangle is a 7.2 barn power-law reference for the data. 
The dashed rectangle is the corresponding participant-scaling 
reference. The sloped dash-dot line is inferred from endpoint 
no of the data and the mean NSD p-p multiplicity. 

Fig. [3] (right panel) shows the same data plotted in the 
power-law format (including the first bin of the analy- 
sis, bounded above by the dotted line at n h - = 5). The 
figure confirms that the power-law format is also appli- 
cable to RHIC data. Quantitative details revealed in the 
linear power-law format are not accessible in the semi- 
log format of the left panel. Two of the five power-law 
parameters can be obtained directly from the data: the 
upper endpoint n^ 4 ~ 4.15 = 297 1 / 4 (mean n h - for 
b = 0) and the fluctuation width (slope) at no. The 
shape about the upper endpoint indicates that fluctu- 
ations are nearly symmetric about hq, in contrast to 
the significant skewness of Hijing (Fig. [5]- left panel). 
Lower endpoint n p and the fluctuation slope at n p are 
not directly accessible due to efficiency and background 
uncertainties and the binning scheme (linear on n h -). 
However, the left edge can be sketched (dotted line at 
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left endpoint) based on the expected p-p h~ yield. The 
corrected n^^ = n NNh + + n NN h - at 130 GeV should 
be ^2.3 in one unit of rapidity [la ]. Thus, we estimate 
np /4 « (0A5n NN , h -) 1/4 = 0.85. 

The solid rectangle represents the power-law reference 
corresponding to <tq = 7.2 barns distributed uniformly on 
the interval between the estimated n p ^ 4 and the observed 
n^ 4 . The dashed rectangle represents the participant- 
scaling reference, where n^ 4 — > {n^) 1 ^ 4 = (^wiV./i- ' 
n P arts)l'*) 1/i = (1.15-191) 1 / 4 = 3.85. Combining the two 
references we can predict the average slope of the data 
distribution between endpoints. Since the total cross sec- 
tion is the same for both data and participant-scaling 
reference the slope is m — — 2/l(n 1 / 4 — [n-o] 1 ^ 4 ) / '( n c/ 4 — 
[n N N,h-} 1/4 ) 2 = -0.15, where A = da/dn 1 ^ ~> 2.4 
barns is the uniform differential cross section for partic- 
ipant scaling, with n\ , riNN,h- and ( n o) 174 evaluated 
above. The dash-dot line drawn with that slope starting 
at {n^N .h- ) 174 forms a trapezoid approximation to the 
data. The fractional change in height is 0.15 x 3.3/2.2 « 
22%, consistent with two-component parameter x = 0.08, 
as shown in App. |XJ 

The comparison with data indicates good agreement 
between the results of [l2j and the power-law description 
with N-N (p-p) constraint. The two equivalent physics 
results from the analysis are the negative slope of the 
power-law distribution and the difference between no val- 
ues for data and the participant-scaling reference (c/. 
App. lC 2| ). In essence, given endpoint n and the NSD N- 
N (p-p) multiplicity the entire minimum-bias distribution 
is known sufficiently well for centrality determination at 
the percent level. This review of a conventional central- 
ity analysis in a power-law context provides some idea of 
the precision possible with the power-law format. 



VI. THE GLAUBER MODEL 

The Glauber model of nucleus-nucleus collisions rep- 
resents the multiple nucleon-nucleon interactions within 
the two-nucleus overlap region in a simply calculable 
form. The nuclear-matter distribution is modeled by 
a Woods-Saxon (W-S) function. The nucleon distribu- 
tion can be modeled as a continuum W-S distribution 
(so-called optical Glauber) or as a random nucleon dis- 
tribution sampled from the W-S density (Monte Carlo 
Glauber) . The Glauber geometry parameters are partic- 
ipant number n par t, N-N binary-collision number nbim 
nucleon mean path length v = 2 nun/n par t and A- A cross 
section a{b). Precise determination of n par t, nu n and v 
in relation to impact parameter b and fractional cross 
section a(b)/ao establishes the centrality dependence of 
particle production and correlations. The following de- 
scriptions of optical and Monte Carlo procedures are de- 
rived in part from [l3l . [Llj . 



A. Optical Glauber 

The optical Glauber relates A-A geometry parameters 
to b through continuous integrals of the nuclear den- 
sity. Normalized function pa^) is a 3D nuclear den- 
sity with a Woods-Saxon radial form. The projection 
onto a plane normal to collision axis z (single-particle 
areal density) is defined by Ta(s) = J dz Pa{s, z). The 
overlap integral (two-particle areal density) for nuclei 
A and B is ?ab(&) = / ds Ta{s)Tb{s — b), an auto- 
correlation distribution for Ta(s) if A — B. n part (b), 
nbin{b) and da(b)/dirb are defined as integrals of Ta 
combined with the appropriate nucleon-nucleon cross sec- 
tion a NN [13]. n bin (b) = J dsAT A {s) a NN BT B {s - b) 
and for A, B ^> 1, da(b) / dub 2 ~ 1 — exp{— nu n (b)} and 
n part (b)/2 ~ / dsAT A {s){l - exp[~a NN BT B (s - b)]}. 
The expression for da(b)/dirb 2 implies that n&i n (6o) ~ 
ln(2) < 1 (and actually goes to zero for large b), whereas 
we expect ntm —* 1 from above for peripheral collisions 
of real nuclei. The same problem arises for n par t, since 
npart/2 — > n bin for large b. 

B. Monte Carlo Glauber 

The Monte Carlo Glauber simulates an ensemble of 
A-B nucleus-nucleus collisions for a distribution of im- 
pact parameters b uniform on b 2 . Each simulated colli- 
sion combines discrete nucleon distributions sampled ran- 
domly from the continuous Woods-Saxon nuclear den- 
sities pa(j), Pb{^)- From the event ensemble n par t{b) 
and n b inib) are sampled as correlated random variables. 
Minimum-bias differential cross-section distributions on 
n P art and nbin are constructed from those data. 




100 200 300 400 500 1000 

n part n bin 



FIG. 4: Left panel: The conventional minimum-bias distri- 
bution on participant number obtained from a Monte Carlo 
Glauber simulation [ill ]. Right panel: The corresponding dis- 
tribution on binary-collision number. 

Fig.[4]shows minimum-bias distributions from a Monte 
Carlo Glauber simulation for 200 GeV Au-Au colli- 
sions [IH plotted in the conventional semi-log format. 
When plotted in a log- log format the distribution on n par t 

closely follows the power-law trend n p ^/ 4 similar to data, 
whereas the distribution on n^m follows the power-law 
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trend n 



-5/6 
bin 



Those trends suggest transformations to 



power-law distributions on n p ^ t and 



C. Power- law Glauber 

In Fig. [5] we plot power-law minimum-bias distribu- 
tions on (npart/2) 1 / 4 and n\(^ from the Monte Carlo 
Glauber data in Fig. [4] The distributions are nearly 
rectangular, and bounded on the right end by endpoints 
n pa rtfi/2 = 191 and n binfl = 1136, with (n part;0 /2) 1/4 = 



3.72 and n 



1/6 
6m. 



3.23. The lower endpoints n par t iP /2, 
nun,p an d binning scheme (dotted lines) are discussed 
in the next section. The distribution on (npart/2) 1 / 4 is 
especially simple: a constant plus a sinusoid with 5% 
relative amplitude described by 



da 



d(n part /2)^ 



1/4 
J part,0 



1/4 

7 

part,p 



(i) 



|l + 0.05 sin ^2.6 (n part /2) 1/4 - 2.3 U 



with centroid 2.3 = {(n part ,o/2) 1/4 + (n part , p /2) 1 / 4 }/2 = 
(191 1 / 4 + 0.5 4 / 4 )/2. That expression is plotted as the 
dashed curve just visible at the left end of the left panel. 




FIG. 5: Power-law minimum-bias distributions on 
participant-pair (left panel) and binary-collision (right 
panel) numbers, obtained from a Glauber Monte Carlo 
simulation [ill ]. The assumed total cross section is 00 = 7 
barns. The upper endpoints (half-maximum points) for 

rip{ 4 t and n\{^ are 3.72 and 3.23 respectively, and the lower 
endpoints are 0.5 1/4 and 0.5 1/6 respectively. 



In App. [B] we discuss several problems in relating 
Glauber parameters to the fractional cross section with 
these differential distributions. A critical issue is how 
the integer spaces should be binned and integrated. For 
the integration method used to obtain dash-dot curves in 
Figs. [6] and [7] the dotted lines in Fig. [5] represent the first 
three bins. The dash-dot curve in the left panel passes 
through the original Monte Carlo data distribution on 
n-part- The solid curve passes through a regrouped distri- 
bution on integer n part /2 consistent with assumptions in 
the participant-scaling hypothesis fApp. |C~2|) . 



The power-law format reveals details at the percent 
level inaccessible with the semi-log format of Fig. |4j The 
observed approximate power-law trend for the minimum- 
bias distribution on n c h is a consequence of the nearly- 
exact power-law trend on n part . Whatever its origins, 
we capitalize on the simplicity of the power-law trend to 
refine centrality measurement and better understand the 
mechanisms of particle, p t and E t production. 



VII. GLAUBER PARAMETERIZATIONS 

Figs. [H] and [7] show running integrals of Monte Carlo 
Glauber data which connect n part /2 and n i n to central- 
ity measured by the fractional cross section in the form 
1 — o-/o~q. For the purpose of centrality determination (to 
~ 2%) the power-law Glauber curves at 200 GeV are well 
represented by simple linear expressions (npart/2) 1 / 4 = 
ct/cto • (n part , p /2) 1 / i + (1 - cr/<7 ) • (n-part,o/2) 1/4 and 



i 1/6 



■ n b(n,p + (1 - 0-/0-0) 



l 1/6 



The dotted 

lines in Figs. [6] and [7] represent those power-law refer- 
ences, with lower endpoints n par t lP /2 and nbi n ,p set equal 
to 3/4 (upper dotted lines) and 1/2 (lower dotted lines). 
The correct endpoint choice (1/2) is justified in App. [Cl 
The dashed and dash-dot curves are running integrals 
of the distributions in Fig. [5] (cf. App. [SJ. Within the 
power-law format we parameterize the dash-dot curves 
simply and precisely as the solid curves. The insets pro- 
vide details of the peripheral regions. 



A. Running-integral definitions 

In Fig.[S]four forms of running integral (two dashed and 
two dash-dot curves) are plotted. The solid curve is the 
parameterization defined in Eq. The dashed curves, 
consistent with the power-law reference with endpoint 
n par t,p/2 = 3/4 (upper dotted line), can be obtained in 
two ways (/$ is the i th element of a minimum-bias differ- 
ential histogram in Fig. [4] with M elements): 

1) Unit bins on continuous n part are centered on inte- 
ger n par t values. The running sum of /, is plotted at up- 
per bin edges on (n par t/2) 1//4 , defining a middle Riemann 
sum on n part . The running sum is F m = Ylj—i fi, and the 
normalized running sum is 1 — ct/cto = F m / Fm- Its m th 
value is plotted at upper bin edge {{n part ^ m + l/2)/2} 4 / 4 
(open squares in Fig. [S]- inset) . The lower endpoint (low- 
est bin edge) on n part /2 is 3/4 {n part , mm = 2), and the 
bin edges continue upward on odd quarters. 

2) Unit bins lie between integers on continuous n par t- 
At the i th integer value of n part the differential cross- 
section entry at that and the proceeding integer n part 
value are averaged, defining an upper Riemann sum. The 
running sum is then F m — J2iL2(fi + /i— 1)/2, plotted at 
upper bin edge (?T. pQr t im /2) 1 / 4 . (open triangles in Fig. [5] 
- inset) The difference between 1) and 2) is mainly that 
points on (n part /2) 1 / 4 are shifted by 1/2 bin. The lower 
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FIG. 6: Participant-pair number (np ar t/2) 1//4 vs fractional 
cross section 1 — cr/aa in the power-law format. The dotted 
lines are power-law references for two choices of lower end- 
point (0.5 1//4 and 0.75 1 / 4 ). The dash-dot curves are running 
integrals of the distributions in Fig. [5] (left panel, solid and 
dash-dot curves) with endpoints n pa rt,p/2 = 1/2. The dashed 
curves represent alternative running integrals with lower end- 
point n part ,p/2 = 3/4. The solid curves are the n par t/2 pa- 
rameterization described in this section. The solid dots are 
from a published Monte Carlo Glauber simulation [l4l ]. 

cndpoint is also 3/4. Either method emulates bin aver- 
ages on rip ar t used to relate n part to u/ctq. 

The dash-dot curves, consistent with the power-law 
reference with endpoint n part .p/2 = 1/2 (lower dotted 

line) can also be obtained in two ways (/* is the i th 
element of the power-law minimum-bias differential his- 
togram derived from /j, with M elements): 

1) Unit bins lie between integers on continuous n part . 
Each entry from the power-law histogram on n pa 4 t is 
multiplied by the bin width on (npart/2) 1 ^ 4 preceding 
it. The running sum is plotted at the (n part /2) 1 ^ 4: value 
of the entry, defining an upper Ricmann sum on n par t- 
The running sum is F m = YJI^fl^ ' [{ n P art.i/2) 1/i - 
(npart^-i/2) 1 / 4 ], and the corresponding normalized run- 
ning sum is plotted at (n pa rt,m/2) 1 ^ 4 (open circles in 
Fig. [5]- inset). The effective lower endpoint (lowest bin 
edge) on n part /2 is 1/2. 

2) Unit bins on continuous n part /2 are centered on 
integer values oin par t/2. The first few bins are illustrated 
by the dotted lines in Fig. [5] (left panel). The M entries of 
the f 1 ^ 4 power-law histogram on n^ 4 rt transformed from 
fi on n par t (dash-dot curve in that panel) are combined 

in pairs to form M/2 entries g^ 4 on integer values of 
n P art/2 (solid curve in that panel), with g 4 J 4 — f 1 ^ 4 + 
/^j+i- The corresponding bin widths on (npart/2) 1 / 4 are 
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FIG. 7: N-N binary collision number n\{^ vs 1 — er/oo. The 
dotted lines are power-law references for two choices of lower 
endpoint (0.5 1/6 and 0.75 1/6 ). The dash-dot curve is a run- 
ning integral of the distribution in Fig. (right panel) with 
lower endpoint n\H = (1/2) 1/6 . The open circles in the in- 
set are the first few points on the dash-dot curve. The solid 
curves are the nn„ parameterization described in this section. 
The solid points are from a Monte Carlo Glauber simulation 
with binned centrality [l4| . 



determined exactly, and the running sum is plotted at 
upper bin edges on n part /2, defining a middle Riemann 

sum on n part /2. The running sum is G m = Y^=i ' 
[(n pal . t , J /2+l/2) 1 /4_( npart ^./ 2 -l/2) 1 /4] i an d G m /G M/2 

is plotted at (n par t , m /2 + 1/2) 1 / 4 (solid dots in Fig.©- 
inset). The effective lower endpoint (lowest bin edge) on 
n par t/2 is again 1/2. This last method complies fully 
with the participant-scaling hypothesis (cf. App. IC 2[) . 
Bin averages on n part /2 should be consistent with these 
results. 



B. Full power-law parameterizations 

The simple linear power-law parameterizations de- 
scribed in the beginning of this section provide good vi- 
sual comparisons and specific expectations for limiting 
cases (extrapolation constraints). For particle production 
studies the accuracy demand on n part /2 increases sub- 
stantially, and the 5% sinusoid should be included in the 
[npart/2) 1 ' /4 parameterization, as discussed in Sec. IIX Al 
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The (ripart/2) 1 / 4 parameterization is then elaborated to 



/ Tlpart \ 

\ 2 ) 



1/4 



G / Tip ar i p \ ^/^ 

a \ 
0.05 



2.6 
0.05 

~2lT 



2 

■ cos ( 2.6 



<j j \ 2 
(n' part /2)^-2.3 



(2.6[(n partiP /2) 1 / 4 -2.3]), 



(2) 



where (n' part /2) 1 / 4 is defined by the linear parameteriza- 
tion (first line), n parU) /2 = 191 at -y/s = 200 GeV and 
n P art, P /2 = 0.45. 

The solid and dash-dot curves in Fig. [5] agree well, de- 
viating from the lower dotted reference line (endpoint at 
n par t/2 = 1/2) by a small curvature concave downward 
due to the sinusoid component. In the peripheral region 
(see inset) the thin solid line (parameterization without 
fluctuations) agrees almost exactly with the thick dash- 
dot line and solid dots from method 2) of the power-law 
integration, which fully reflects the participant-scaling 
hypothesis. The full parameterization (thicker solid line 
in inset) includes an accommodation for n c h fluctuations 
[Eq. (j4j] such that n part /2 — > 1 for peripheral collisions 
(c/. App.EJ. 

The solid and dash-dot curves in Fig. [7] also agree 
well. The definition of running integration is simpler for 
ribin, since the middle Riemann sum on n&i n is consistent 
with binary-collision scaling and bin averaging. Because 
the da/dnl(® distribution in Fig. [5] (right panel) begins 
well below the mean value and has a significant positive 
slope, the corresponding running integral in Fig. [7] (dash- 
dot curve) has a significant curvature concave downward 
which we accommodate by a modification of the linear 
power-law reference, 



i 1/6 

bin 



, I \0.965 
(<7/<7 ) " 



liL + (1 - t/^o) 



1.035 1/6 

' "f>m,0- 



(3) 



The exponents on the cross-section factors add the neces- 
sary curvature to the parameterization, as shown by the 
close agreement between solid and dash-dot curves. 

The final parameterizations in Figs. [5] and [7] agree to 
~ 1% with the power-law integrals of the Glauber Monte 
Carlo data, except for the peripheral region where effects 
of multiplicity fluctuations arc modeled in the full param- 
eterization. To obtain the asymptotic approach to unity 
required for n par t/2 and nun in peripheral A- A and N-N 
collisions we use 



(1 



a\l/a 



(4) 



,1/6 



with a = 8 for x = {npart/2) 1 ! 4 and a — 15 for x = nl 
The choices for a depend on observed N-N multiplicity 
fluctuations as discussed in App. [Cl 

Also plotted in Figs. [5] and [7] are results from an 
independent Monte Carlo Glauber analysis (nine solid 
points) based on bin averages on n part and nu n • For 
{npart/2) 1 ^ in Fig. [5] the results agree with the power- 
law parameterizations for the most central collisions but 



significantly disagree for more-peripheral collisions. The 
points are consistent with the linear parameterization of 
(npart/2) 1 ^ 4 with endpoint n par t tP /2 = 3/4 (upper dot- 
ted line) . The dashed curves which pas through the solid 
points represent alternative running-integral definitions 
on n part described above, with endpoints at (3/4) 1 / 4 . We 
conclude that the incorrect endpoint is a consequence of 
bin averaging on n part rather than n par t/2. 

The rc^f data (solid points) from [l4| in Fig.[7Jseem to 
be consistent with endpoint n i n ,p = 3/4 (upper dotted 
line). However, although the points are slightly higher 
than the solid and dash-dot curves, both data and curves 
are displaced from the lower dotted line with endpoint 
1/2 by a curvature resulting from the structure of the 
differential cross section of Fig. [5] (right panel). Thus, 
the ribi n data points are probably consistent with the 
correct endpoint 1/2. 



C. Power- law mean path length v 

We have obtained precise parameterizations for 
{npart/2) 1 ' 14 and vs fractional cross section. We now 
define participant path-length estimator v. The path- 
length concept originated with h- A experiments [la ] , for 
which v was defined in terms of a hadron interaction 
length in nucleus A depending on the hadron-nuclcon 
cross section and center-of-mass energy. The correspond- 
ing definition based on Monte Carlo Glauber parameters 
is v = 2nbi n /n par t- Using the full power-law parameteri- 
zations from Eqs. @ and we have 



v(a/a ) 



1/6 



2n bi . 



npart ([n part /2]i/4) 



(5) 



In Figs. [6] and [Jj the n par t/2 and n i n parameterizations 
(solid curves) which accommodate n c /j fluctuations are 
constrained by Eq. (g]) so that n part /2 and n bin > 1 as 
1 — ct/(7q 0. Those transitions must be coordinated so 
that v — * 1 smoothly as well. 

Fig. [8] (left panel) shows v(o /<jq) from Eq. ([5]) as the 
solid curve. Maximum value vq ~ 6 corresponds to b = 
and &nn f° r 200 GeV N-N collisions (n oin cx <Jnn f° r 
central collisions), vnn ~ 1-25 corresponds to the cen- 
trality (~ 0.06) at which n part /2 ~ n^m ~ 1 (cf. solid 
curves, right panels). The dotted curve (barely visible 
near the center) shows the result when the sinusoid is 
omitted from the n par t/2 parameterization. The dash- 
dot curve (ripart/2) 1 / 3 , a fair approximation to v except 
in the most peripheral region, is used to approximate v 
in Eq. (|A3[) . The dashed curve which underlies the solid 
curve results from changing both lower endpoint values 
for Eq. ([5]) from 1/2 to 3/4. The changes produce no 
significant change in showing that v provides strong 
common-mode reduction of sensitivity to binning and av- 
eraging schemes. 
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FIG. 8: Left panel: Participant mean path length v vs frac- 
tional cross section. The solid curve is derived from the 
parameterizations of the running integrals in Figs. [6] and [7] 
(solid curves), and the dotted curve is the same but with- 
out the n par t sinusoid. The dash-dot curve is approximation 
(n V art/2) 1 ^ 3 . The points represent two binned Glauber re- 
sults [lj] . Right panels: Detailed views of Glauber power-law 
trends for peripheral collisions. The solid lines with endpoints 
at (1/2) 1/4 and (1/2) 1/6 are derived from the dash-dot run- 
ning integrals in Figs. [6] and [7] The points are from opti- 
cal (open) and Monte Carlo (solid) Glauber simulations [14| . 
The dashed lines extrapolate the solid points to endpoints at 
(3/4) 1 / 4 and (3/4) 1 / 6 . 

The extrapolation endpoints of the linear n part /2 and 
nbin power-law trends are both ~ 1/2. However, the 
physical limit for each parameter is 1. Resolution of that 
conflict is addressed in App. [C] It implies however that 
v = 1 does not represent the point on v corresponding 
to a single N-N collision. To locate that point we make 
the following argument. The fractional cross section cor- 
responding to rip ar t/2 = 1 in the linear power-law case is 
given by 

I vqJnn (v<,o/2) 1/4 -(v<,p/ 2 ) 1/4 

With n par t,p/2 ~ 1/2 and n pa rt,o/2 = 191 we obtain 

< 1 — — > ~ 0.06. vnn ~ 1-25 is therefore the loca- 
L a o j NN 

tion on v of the centroid of the single N-N multiplicity 
distribution, as noted in Fig. [5] and subsequent plots. 

Participant mean path length v is the ideal centrality 
measure for A-A collisions, providing precise visual tests 
of N-N linear superposition relative to a two-component 
combination of participant and binary-collision scaling. 
The per-participant yield of n c h, Pt or E t plotted vs v, for 
a simple combination of participant and binary-collision 
scaling, should exhibit a linear increase with v relative 
to a constant background [l7| (cf. App. [5]). Deviations 
from linearity would reveal A-A medium effects, of cen- 
tral importance to RHIC physics. 

In contrast, the conventional centrality measure 
n par t ~ 2t/ 3 is a very nonlinear centrality measure, and 

1/3 

binary-collision trends emerge as n part curves. As il- 
lustrated in Sec. IX Bl the nonlinearity of n part relative 
to the fractional cross section compresses much of the 



cross section into a small region near zero. The lower 
50% of the cross section occupies less than 15% of the 
n P art range. Central collisions dominate the n par t plot- 
ting format. Thus, the opportunity to discern subtle but 
important deviations from binary-collision scaling start- 
ing from N-N collisions (e.g., jet quenching and other 
medium effects) is abandoned. 

D. Comparison of optical, Monte Carlo and 
power-law Glaubers 

Fig. [8] (left panel) also compares optical and Monte 
Carlo Glaubers with the power-law parameterization for 
v in Eq. (J5])- The open circles and solid points are 
bin mean values from optical and Monte Carlo Glauber 
simulations respectively, presented in Tables II and III 
of [HI]. The agreement is notable in view of the sig- 
nificant discrepancies between Glauber implementations 
in n P art and ribin in the right panels. The solid curve 
from the power-law parameterization also agrees well 

with v = 2(nl(n) 6 / (rip^t) 4 , since those quantities are 
nearly linearly related to the fractional cross section. As 
noted above, v provides excellent common-mode reduc- 
tion of systematic errors. 

In Fig. [5] (right panels) the lower half of the central- 
ity range is shown to increase sensitivity in the critical 
peripheral region. The optical Glauber quantities go to 
zero asymptotically by definition. The solid curves are 
the linear parameterizations (without fluctuations) de- 
rived from the running integrals of minimum-bias dis- 
tributions in Fig. [5] and El with endpoint (1/2) 1 / 4 . The 
dashed curves which agree with the solid points from [l4[ 
have endpoints (3/4) 1 / 4 . The endpoint difference and the 
preferred endpoint choice are discussed in App. [Cj The 
difference between dashed and solid curves propagates to 
a ~ 50% error in n part /2 for peripheral collisions. 

VIII. CENTRALITY DETERMINATION 

We now return to the problem of relating measured 
data to A-A collision geometry. The relation between 
measured n c h and fractional cross section ct/cto is defined 
by running integration of the minimum-bias data distri- 
bution. The multiplicity is in turn related to the Glauber 
geometry parameters through the fractional cross sec- 
tion. The power-law format, running integrals and ex- 
trapolation constraints for N-N collisions greatly reduce 
systematic uncertainties in collision geometry, especially 
for peripheral collisions. 

A. Conventional Centrality 

In a conventional centrality determination the raw 
minimum-bias distribution dN evt /dn c h is corrected for 
event-trigger inefficiency (typically a small effect for all 
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n c h)j vertex-reconstruction inefficiency (possibly large for 
peripheral collisions and small n c h), tracking (particle- 
detection) inefficiencies and backgrounds {e.g., beam-gas 
collisions and photo-nuclear excitations). Tracking in- 
efficiencies depending on n c h distort the minimum-bias 
distribution (e.g., change the average slope of the power- 
law format from negative to positive), and may produce 
systematic errors in the centrality determination. The 
corrected and normalized distribution da(n c h)/dn c h or 
da(n c h)/dn 1 J h 4 should integrate to the total cross section 
cr corresponding to the trigger definition. 

As an example, we consider an analysis of 130 GeV 
RHIC data [14[ ■ We work backward from centrality bin- 
edge definitions to the running integral to the differential 
power-law distribution. Ten centrality classes were de- 
fined on total multiplicity n c h = n h + + n h - detected in 
one unit of pseudorapidity. The raw minimum-bias dis- 
tribution dN evt /dn c h was corrected for significant trigger 
and vertex-finding inefficiencies below n c h = 50 as fol- 
lows. The raw data were scaled to agree with Hijing in 
the trusted n c h interval [50,100] (compare to the [5,25] 
interval in [l^]) where the collision dynamics were said 
to be dominated by A-A geometry and well-described by 
the Hijing model. 

The trigger/vertex efficiency below n c h — 50 was de- 
termined by the ratio of data to Hijing minimum-bias 
distributions, and the data were corrected by that ratio, 
giving a reported overall trigger/ vertex efficiency of 94%, 
with a 60% efficiency below n c h = 5. The corrected dis- 
tribution was normalized to the 6.9 barns total cross sec- 
tion do at 130 GeV estimated from Hijing and partitioned 
into ten bins according to the integrated fractional cross 
section. The total of beam-gas and photonuclear contri- 
butions to the most-peripheral 20% bin was estimated 



to be 30% according to [14|, and that bin was therefore 



excluded from further analysis. 

In Fig. [5] (left panel) the points represent reported 
centrality bin edges on uncorrected n c h in \rj\ < 0.5 
and p t > 0.1 GeV/c from Table I of 0. We infer 
upper endpoint no = 488 for uncorrected data by ex- 
trapolating the bin-edge positions with a parameteriza- 
tion (dotted curve). We estimate the lower endpoint as 
Tip — rip ar t , P /2 • n NN — 0.45 x 2.3 = 1.04 for corrected 
n c h based on the expected 130 GeV p-p yield in a pseu- 
dorapidity acceptance of one unit [lg]. The TPC track- 
ing efficiency was reported to be 80%, corresponding to 
n p = 0.83 for uncorrected data. The dashed line repre- 
sents the participant-scaling reference for corrected data 
and the solid curve estimates the power-law reference for 
corrected data. 

In Fig. [9] (right panel) the dotted lines and curve recon- 
struct the uncorrected differential cross section by differ- 
entiating the parameterization of the data points (dotted 
curve) in the left panel. The dashed rectangle represents 
the participant-scaling reference for corrected data, based 
on tinn = 2.3 and (To = 6.9 barns, and the solid rectan- 
gle represents an estimate of the power-law reference for 
corrected data with no = 1.4 n par t,o/2 • unn = 615. The 
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FIG. 9: Left panel: The points are bin edges from a con- 
ventional centrality definition for uncorrected particle multi- 
plicities n c h [lj. The dotted curve is a parameterization of 
the points. The dashed curve is the participant-scaling refer- 
ence for corrected n c h- The solid line (power-law reference) 
is an estimate of the corrected data trend for 130 GeV. Right 
panel: The corresponding power-law form of the differential 
cross section inferred from data (dotted curve and lines), the 
participant-scaling reference (dashed rectangle) and the esti- 
mated power-law reference for corrected data (solid rectan- 
gle). 



dash-dot line estimates the expected trend of the cor- 
rected differential cross section and, given the 2x factor 
for h~ — > n^, compares well with Fig. [3] (right panel) for 
the h~ yield from the same 130 GeV data [121 ]. 

The uncorrected distribution represented by the dotted 
lines and curve from [1 41) differs markedly from the cor- 
rected results from [12| represented by the dash-dot line. 
The reason is a strong n c ^ dependence of the tracking 
efficiency. The estimated 80% TPC tracking efficiency is 
not uniform on n c /j. A TPC tracking efficiency typically 
decreases significantly with increasing n c h due to losses 
by cluster/track merging for larger track density. The 
result is the difference between uncorrected data (points 
and dotted curve) and estimated corrected data (solid 
line) in the left panel. The inefficiency variation reverses 
the curvature of the running integral in the left panel and 
changes the sign of the slope at mid-centrality in the right 
panel (cf. App. |A"|) . 



B. Power-law Centrality 

Power-law centrality determination employs a running 
integral of the power-law differential cross-section distri- 
bution, removing the restriction to fixed centrality bins 
and a special n c h definition. The general features of 
the running integral are represented by the cartoon in 
Fig. [TO] (left panel) . The participant-scaling reference is 

the dashed line running from n\j 4 ~ (n pQrtjP /2 • njvAr) 1 ^ 4 
to (n^) 1 / 4 = (n partj o/2 • njvAr) 1 / 4 on the vertical n*^ 4 
scale and from to 1 on fractional cross section 1 — (t/cto- 
The power-law reference is the thin solid line running 
from n\j 4 to n\j 4 . Those lines correspond to dashed and 
solid rectangles respectively in previous differential cross- 
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section plots. 

Endpoints n p , no and the slope at the midpoint of 
the differential cross section determine the shape of the 
running integral (solid curve) . The differential endpoints 
determine the mean slope of the running integral. The 
differential slope determines the curvature of the running 
integral at its midpoint. The slopes at the ends of the 
differential distribution (from fluctuations) determine the 
curved segments at the ends of the running integral. The 
fluctuation shape near n p also depends on the skewness 
of the N-N multiplicity distribution, which causes a shift 
of n p slightly below n^^/2. 




FIG. 10: Left panel: A cartoon of a running integral of 
the power-law differential cross section showing the principal 
features. Right panel: Running integrals of the differential 
cross sections in Fig. [2] (left panel) for quench-off (dash-dot) 
and quench-on (thicker solid) Hijing. The dashed line is the 
participant-scaling reference. The thinner solid line is the 
power-law reference. 



Fig. [10] (right panel) shows running integrals of the 
Hijing quench-on (solid) and quench-off (dash-dot) dif- 
ferential power-law distributions in Fig. [5] (left panel). 
Running integration substantially reduces statistical 
noise 18]. Dashed line ■ {n part /2) 1 / 4: represents 

participant scaling (particle production proportional to 
participant-pair number). Deviations from the power- 
law reference (thin solid line) should consist mainly of 
the small curvature in the central region corresponding to 
the slope of the differential distribution. Excursions near 
the endpoints correspond to multiplicity fluctuations. 

The curvatures for the Hijing data are concave upward, 
corresponding to the negative slope of the differential dis- 
tributions in Fig. [5] (left panel). Given that n p ~ unn/2 
(exact for a symmetric N-N multiplicity distribution), 
knowledge of the NSD p-p multiplicity to ~10% means 
that the centrality accuracy for peripheral A-A collisions 
is ~ 1%. For Hijing, the relation is n p ~ QA5unn = 2.25 
due to the skewness of the N-N NBD distribution (c/. 
App. p. In Fig. [10] (right panel) n p = 0.45 • 5 = 2.25 re- 
sults in excellent agreement between Hijing data (curves) 
and the participant-scaling reference (dashed line) for pe- 
ripheral collisions. 



IX. PARTICLE, Pt AND E t PRODUCTION 

One purpose of A-A centrality determination is to es- 
tablish how various quantities are 'produced' in the fi- 
nal state from a multitude of initial N-N (or parton- 
parton) collisions, and how (or if) the production mech- 
anisms change with A-A collision geometry. Produc- 
tion of final-state hadrons, transverse momentum p t and 
transverse energy E t is determined by initial-state par- 
ton scattering, parton dissipation in the bulk medium, 
the dynamics of the bulk medium and the hadronization 
process (scattered-parton and bulk-medium fragmenta- 
tion). Change can be defined relative to a linear superpo- 
sition hypothesis: a linear combination of products from 
the number of N-N collisions predicted by combining the 
Glauber and two-component models. 

Tick, Pt and E t per participant pair should exhibit a 
combination of participant scaling (independent of v) and 
binary-collision scaling (proportional to v) , plus medium 
effects which may be substantial. The differential proce- 
dure consists of two parts: 1) determine the ratio of a pro- 
duced quantity to the number of N-N participant pairs, 
and 2) plot the variation of that ratio vs path length is, 
the number of N-N collisions per participant pair. That 
procedure can be illustrated with the Hijing Monte Carlo 
and RHIC data. In contrast to particle number produc- 
tion, the pt and E t minimum-bias distributions involve 
continuous variables which require some differences in in- 
tegration technique. In the case of E t production we 
compare conventional and power-law methods applied to 
RHIC data. 



A. Particle production 

The two-component model of nuclear collisions [13, Gil 
compares the fraction of particle production in A-A colli- 
sions due to participant scaling n c h cx n part /2 (soft com- 
ponent) and the fraction due to binary-collision scaling 
rich of n bin (hard component). That decomposition may 
separate contributions from initial-state parton scatter- 
ing and fragmentation, which should scale as the latter, 
from bulk-medium hadronization which may scale as the 
former. By comparing deviations from participant scal- 
ing with a linear binary-collision reference on participant 
path-length is, modifications to parton scattering and 
fragmentation by the QCD medium (e.g., jet quenching) 
may be analyzed. 

Measurement of particle production in the form 
^■/n P arfTich is a stringent test of centrality determination, 
since the relation of n c h to n par t/2 is sensitive to small 
relative errors in a/ao inferred separately in the data and 
Glauber contexts. Conventional centrality methods with- 
out extrapolation constraints entail uncertainties for pe- 
ripheral collisions large enough that particle-production 
studies with conventional methods have not been at- 
tempted for the 20% most peripheral collisions. The 
power-law method opens that region to precise study. 
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In the simplest power-law method n*{^ and n^l® vs 
<J /<Jo are represented by linear parameterizations suffi- 
ciently accurate for centrality determination. However, 



the sinusoid on dajdnj 4 rt is 0.05 sin{2.6( 



i 1/4 

part 



2.3)} 

(relative to the mean power-law cross section). The re- 
sulting relative deviation of n p ^. t vs ^(nl/^) / ao from 
a power-law trend is 0.014 at the midpoint, an error 
of 1.5% in centrality determination. However, particle- 
production studies involve n part /2 and not (npart/2) 1 / 4 . 
The sinusoid contributes a deviation in n part /2 of 6% 
from the power-law reference, sufficient to require includ- 
ing the sinusoid in particle-production studies. There- 
fore, the sinusoid in Fig. [5] (left panel) should be included 

in the parameterization of rhj 4 rt vs c/frj. 

Fig. [IT] (left panel) is a cartoon of 2/n par t ■ n c h vs 
mean participant path length v that combines nH^ - 



dot) and quench-on (solid) Hijing. The dotted curve 
shows the result of omitting the sinusoid from the n par t 
parameterization — a 6% systematic error for particle pro- 
duction. Parameter a in the Glauber parameterization 
which represents fluctuations affects the curves only near 
v = 2. For 1/a — > the slopes of the curves become neg- 
ative near v — 2 because n par t/2 without fluctuations 
decreases too rapidly with decreasing v. The adopted 
value a = 8 for n part /2 was the largest value for which 
the slopes of the curves were non-negative near v = 2. 

The hatched region (mean value and error band) esti- 
mates participant scaling from dn c hjdr\ ~ 2.5 for NSD 
p-p collisions {v w 1.25) at 200 GeV integrated over 
two units of pseudorapidity, consistent with SPPS p-p 
data [2(J. The width of the error band represents the 
systematic uncertainty propagated from the 10% uncer- 
tainty in n partiP j2 = 0.45 ± 0.05 which produces a 5% 



1 - a/a obtained from a data minimum-bias distribution maximum uncertainty in 2 /n part dn ch /dr) confined to the 



with n~J 4 rt and v vs 1 



<t/<7o obtained from a Glauber 
Monte Carlo. The lower band represents data uncor- 
rected for tracking inefficiency (including n c h dependence 
which produces the curvature), and the upper band rep- 
resents corrected data, unn is the N-N (^NSD p-p) 
multiplicity in the acceptance derived from separate ex- 
periments, and e is the tracking efficiency for peripheral 
collisions. The horizontal line extending to vo (b = 0) 
is the participant-scaling reference. The curves represent 
deviations from power-law trends due to fluctuations. 




•5 8 




I 1 | I l 1 | I 1 ! 




a 7.5 




! 1400/191 j 




1 7 








1 6.5 




j 1210/191 / , 




6 




'< quench on / 




5.5 




I ..-^^T- ''quench off 




5 




\pi\\\\\\ii 




4.5 




/ ">NN 




4 




AN-N A-A, b=0 




3.5 








3 









FIG. 11: Left panel: A cartoon of particle production per 
participant pair vs mean participant path-length v for uncor- 
rected (lower) and corrected (upper) data. The effects of mul- 
tiplicity fluctuations are illustrated at the ends. Right panel: 
Particle production curves for Hijing quench-off (dashed) and 
quench-on (solid) collisions. The dotted curve illustrates the 
consequence of omitting the sinusoid in the n part /2 param- 
eterization (linear approximation). The hatched region rep- 
resents the N-N (NSD p-p) multiplicity in the acceptance at 
200 GeV (njvjv ~ 5 for Pythia). The horizontal solid lines 
represent 2no/n par t,o- 

Fig. [11] (right panel) shows the result of a particle- 
production study of the Hijing Monte Carlo, with 
n par t(b)/2 vs a(b) derived from the Glauber parameteri- 
zation with sinusoid (solid curve) in Fig. [6] and n c h vs 
<j(n c h) derived from the solid and dash-dot curves in 
Fig. [TU] (right panel). We plot 2/n par tAr]dn c h/dr), with 
Arj = 2. The curves correspond to quench-off (dash- 



region around v — 2. The curves above v = 2.5 are in- 
sensitive to that endpoint uncertainty. 

Interpretation of the particle-production evolution 
from Fig. [11] (right panel) is straightforward. Below 
v ~ 2 there is participant scaling and fluctuations. Above 
v ~ 2 the quench-off trend increases with path length, 
presumably reflecting increased hadron production from 
(mini-)jets produced in multiple N-N collisions. The 
quench-on trend shows additional increase above v ~ 2.5. 
Presumably, parton energy loss (quenching proportional 
to final-state path length, also oc v) is converted to addi- 
tional hadron production. Unprecedented access to such 
a detailed picture of particle production illustrates the 
importance of the power-law centrality method. 



B. Transverse momentum p t production 

We next study the power-law distribution and p t pro- 
duction. To obtain the distributions in Fig. [2] (left 
panel) the conventional minimum-bias distribution was 
first formed on integer multiplicity n c h, then rebinned 
onto nj^ . For continuous variable pt, in contrast, p]^ 4 
was binned into 50 equal bins, providing adequate res- 
olution for end-point structure. Event numbers were 
accumulated directly into the bins on p]^ 4 (not p t ). 



Fig. [T^l (left panel) shows da/dp]^ 4 vs Pt^ ■ The points 
on the left edge illustrate the uniform bin spacing and 
edge resolution. The upper half-maximum point pto es- 
timates the correspondence on total pt of b — and 
n par t, o/2, and the lower half- maximum point pt p esti- 
mates n parttP /2 ■ pt,NN, roughly half the total p t in the 
acceptance for N-N collisions. The endpoints for Hijing 
are p tp ~ 0.86 - 0.45 p tAW GeV/c and p t0 - 580 (660) 
GeV/c for quench-off (quench-on) events. 

In Fig. [TJ] (right panel) we plot the total per- 
participant-pair quantity 

2/n part A7jdp t /dr] = 2/n part Ar] dn ch /di] ■ (p t ) (7) 
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FIG. 12: Left panel: The power-law form of the minimum- 
bias distribution on total pt within \n\ < 1 for quench-on 
(solid) and quench-off (dashed) Hijing. The points at the 
left end illustrate the binning scheme uniform on p x J A . Right 
panel: The total pt per-participant pair within the acceptance 
vs mean participant path length v for the two Hijing classes. 
The lower solid and dash-dot curves represent twice the per- 
final- state-particle mean p t 2 (p t ) for the two cases. 



within acceptance A77 = 2. The labeled horizontal solid 
lines represent 2ptfi/n par tfi for quench-on and quench- 
off events. The hatched region represents pt,NN ~ 1-9 
GeV/c - n NN -QM GeV/c. The left endpoint p t , p ~ 0.86 
GeV/c in the left panel is about 0.45 x the N-N value in 
the right panel, as expected. The slope for quench-on 
Hijing corresponds to x Pt ~ 0.16, compared to x Hch ~ 
0.08. That difference is qualitatively consistent with the 
trend of increasing (p t ) discussed below: p t increases with 
centrality faster than n c h, and we expect x Pt > i„ c(i . 

The (overlapping) solid and dash-dot curves below the 
hatched region show the ratio of p t production in Fig.[T2l 
and multiplicity n c h production in Fig. 111! Quench-off 
and quench-on events exhibit the same (p t ) variation, 
roughly consistent with trends for RHIC data at 130 and 
200 GeV [H and 200 GeV p-p results (y ~ 1.25) 0. 

Particle and p t production systematics from Hijing col- 
lisions are closely related to (pt) fluctuations and re- 
lated pt correlations from that model representing minijet 
structure [24[. The variation of Hijing per-particle (pt) 
fluctuations with centrality was found to be small [23[ , in 
disagreement with RHIC (pt) fluctuation measurements 
reported in [l[. The weak centrality dependence of Hijing 
minijet-related pt angular correlations is also very 
different from corresponding RHIC data [4|. Since those 
analyses are based on per-particle fluctuation and corre- 
lation measures we conclude that quench-off Hijing scales 
minijet production and total particle production almost 
identically. Hijing quench-off is simply a linear superpo- 
sition of Pythia N-N collisions whose number follows a 
combination of participant and binary-collision scaling. 
In contrast, fluctuation and correlation analysis of RHIC 
data reveals that Au-Au collisions exhibit strong devia- 
tions from such linear superposition. 



C. Transverse energy E t production 

Fig. [13] (upper-left panel) shows a semi-log minimum- 
bias distribution on transverse energy E t measured with 
an electromagnetic calorimeter (EMCal) patch [25] . The 
data were corrected for trigger inefficiencies. The patch 
acceptance was 1 /6 of 27r azimuth over one unit of pseu- 
dorapidity [0,1]. The E t axis was uniformly binned (bin 
width 2.34 GeV) except for the lowest bin (width 1.56 
GeV). Bin sums were converted to densities dN evti i/dE t 
for this analysis, dividing the N evt ^ by bin widths 5E t i 
to obtain the density distribution in the upper-left panel. 

The E t minimum-bias distribution plotted in the 
power-law format in the upper-right panel was obtained 
by multiplying the upper-left distribution bin- wise by Ja- 
cobian 41?^ and normalizing to (To = 7 barns (the de- 
fault value assumed for this study). Because the binning 
is (mostly) uniform on E t the distribution suffers from 
sparse sampling at the low--Et end critical for central- 
ity and E t production studies, and excessive sampling at 
the upper end. Uniform binning on E\^ is preferable, as 
illustrated by the p t analysis in Fig. [T^] (left panel). 

The E t distribution is almost exactly power-law in 
form over the observed E t interval, consistent with par- 
ticipant scaling. The upper endpoint is E^ 4 = 3.24 
= (110 GeV) 1//4 . The lower endpoint must be esti- 
mated. The position of the first data point (roughly the 
half- maximum) at E] /A = 0.78 (GeV) 1/4 establishes the 
rough estimate E tp = 0.37 GeV. However, a lower end- 
point derived from the observed upper endpoint and as- 
suming pure participant scaling is (0.45 • 110/191) 1 / 4 = 
0.26 1 / 4 = 0.71 (GeV) 1 / 4 , with E t , NN = 0.58 GeV also 
inferred. 

Given a common upper endpoint and using the two 
estimates the dashed rectangle in Fig. [T3] (upper-right 
panel) is the participant-scaling reference, and the solid 
rectangle is the power-law reference. The lower endpoint 
should not be higher than that derived from the up- 
per endpoint and a participant-scaling assumption. The 
power-law reference obtained from the lowest data point 
(solid rectangle) is obviously higher than the mean value 
of the data. Since the position of the single data point 
on the lower edge is not a precise estimate we adopt 0.26 
GeV (and the dashed participant-scaling hypothesis) as 
the best model for the data distribution. 

The lower-left panel shows the running integral of the 
minimum-bias data distribution in the upper-right panel. 
The half-maximum points inferred from the upper-right 
panel are represented by the horizontal lines. The diago- 
nal dotted line {E ti NN ■ ^part/2) 1 / 4 with lower endpoint 
at E^ 4 = (0.45 £; t . at at) 1 / 4 is the participant-scaling ref- 
erence. The solid curve and small points plotted at bin 
edges on E^ 4 represent the recommended running inte- 
gral described in App. |Bj The excellent agreement with 
the participant-scaling reference (dotted line) is appar- 
ent. The open triangles are from [25J and also agree well 
with the participant-scaling reference. 
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FIG. 13: Upper-left panel: The conventional semi-log form 
of the minimum-bias density distribution on total Et in an 
EMCal patch (see text). Upper-right panel: The power-law 
form of the normalized minimum-bias density distribution 
on Et for data from 25:]. Lower-left panel: Running inte- 
gral of the distribution in the upper-right panel (solid curve 
and small points) compared to participant-scaling reference 
(n.parf/2) 1 / 4 ■ E t ,NN (dotted line). Lower-right panel: The 
density on r\ of Et per participant pair. The solid and dashed 
curves are derived from the running integral in the lower-left 
panel with n par t /2 parameterizations corresponding to inter- 
cepts n par t,p/2 — 1/2 and 3/4 respectively. The lower solid 
curves compare Et per charge hadron with p t per charged 
hadron. The open triangles in the lower panels are from [2q . 



The lower-right panel shows dEt/dr] per participant 
pair. As noted, a factor 6x was applied to the mea- 
sured EMCal patch E t values to obtain dEt/dr] values 
according to the definition of the patch acceptance. The 
solid curve and solid dots represent the recommended 
integration scheme and N part /2 parameterization from 
the present study. Those results are also consistent with 
participant scaling (hatched region), as expected from 
the uniformity of the power-law minimum-bias distribu- 
tion. Assuming E t .o = E^ , the constant value for all v 
is dEt/dr] = 6 ■ E t ,NN = 3^46 GeV. 

The open triangles are data from [25| ■ The systematic 
error bars in that paper (up to 25% or ~ 0.6 GeV for the 
most peripheral points) are omitted in this comparison, 
since the same underlying minimum-bias data are used to 
compare different integration schemes and Glauber pa- 
rameters. The points from [251 ] (upper triangles) devi- 
ate significantly and systematically from the participant- 
scaling reference (hatched region) and the recommended 
running integral (upper solid curve). The dashed curve 
and open circles represent the correct E^ 4 running inte- 



gral from the lower-left panel and the incorrect N part /2 
parameterization with lower endpoint at 3/4. That com- 
bination describes the points from [25| well, suggesting 
that the incorrect n par t/2 definition with lower endpoint 
at 3/4 contributed to the systematic error in [25j |. The 
systematic uncertainty represented by the hatched region 
is 5%. whereas the uncorrected error in the triangle data 
points approaches 50% at the N-N limit. This example il- 
lustrates the importance of the extrapolation constraints 
available within the power-law context. 

E t and pt centrality trends for data and Hijing are 
said to agree in 25j. However, comparing the consis- 
tent analyses in Fig. [T3] (right panel) and Fig. [T3] (lower- 
right panel) we find that the charged-particle pt and to- 
tal E t (including neutrals) production trends are very 
different. The total E t (hadronic plus electromagnetic 
transverse energy) per participant pair is independent 
of centrality (indicative of participant scaling), whereas 
the number of charged hadrons per participant pair in- 
creases by 40% from peripheral to central collisions at 
200 GeV. Therefore, the total E t per charged hadron 
(e t ) defined in (25| must decrease, whereas the conclu- 
sion in 25j is that (et) increases substantially, similar 
to Hijing n c h and pt production. Solid curves for (et) 
(RHIC data) and (pt)c (Hijing) are compared in the 
lower part of the bottom-right panel. The multiplicity 
trend we used to obtain consistency with [25j at v ~ 6 
is drich/di] — 2.95 (1 + 0.08(f — 1)), whereas we expect 
the prefactor to be 2.5 at 200 GeV. It is notable that the 
solid curves from the present study tend to converge. 

That (et)-(pt)c comparison is still not completely ap- 
propriate because the measured E t represents all par- 
ticles including neutrals, whereas ( vt) represents only 
charge particles. The (et) defined in [25j in terms of n c h 
should instead be defined as (et) = dEt/dr] / dn to t/dr]. 
If we adopt a factor 2 /3 correction from n c h to ntot , as- 
suming hadrons are dominated by pions, we obtain the 
dash-dot curve in the lower-right panel of Fig. [131 and 
the convergence of (et) and (pt)c with increasing central- 
ity toward 0.5 GeV is more apparent. The convergence 
is somewhat artificial (it could be better or worse) be- 
cause Hijing does not necessarily model (pt) from data 
correctly. 

In App. [A] we show algebraically that the four panels 
of Fig. [T5]are redundant. The main experimental result 
of the analysis in [2|| is that E t production in Au-Au 
collisions at 200 GeV and mid-rapidity follows partici- 
pant scaling (xE t ~ 0) within statistical errors. Some 
of the conclusions reported in |25| are in error because 
the Glauber n part /2 trend from [Tj] used to calculate 
2ln par t dEt/d-q is incorrect. The discrepancy is apparent 
in the right panels of Fig. [TBI The correct n part /2 param- 
eterization from the present analysis restores consistency, 
and the physics conclusions change significantly. 
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X. CENTRALITY ERRORS 

We summarize the systematic error sources for two 
centrality methods. We distinguish error sources for the 
Glauber model and for data, and for the conventional 
and power-law methods. The main technical issue is the 
correct relation of data and Glauber parameters to the 
fractional cross section u/cfq. Uncertainties in the relative 
centralities of Glauber model and data, most relevant for 
particle and Pt/E t production, may be significantly re- 
duced by use of extrapolation constraints. 



A. Total cross section and error 

The total cross section for Au-Au collisions is esti- 
mated as follows. The Woods-Saxon matter distribution 
for Au has r = 6.5 fm (including an estimated 0.1 fm 
from the neutron skin) and diffuseness a = 0.5 fm [13| . 
implying a nominal nuclear edge at ro + 2a = 7.5 fm. The 
nucleus-nucleus cross section is then 7r(2 • 7.5 fm) 2 /100 = 
7.05 barns. An uncertainty of 0.1 fm in the edge radius 
results in an uncertainty of 0.2 barns in the total cross 
section. That ±3% range encompasses all Glauber esti- 
mates for RHIC Au-Au collisions at 130 and 200 GeV. 
For comparisons with published data we use stated total 
cross sections. For internal comparisons and illustration 
we adopt the default value <tq — 7 barns. 



B. Glauber- model errors 

We relate Glauber geometry parameters to data 
through their fractional cross-section dependencies. The 
error propagation from fractional cross section to 
Glauber parameters can be estimated from the linear 
power-law parameterizations. For n part , Sn par t/n par t = 
[4(3.72 — 1)](5ct/cto, and the coefficient in square brack- 
ets is 10.9. For nun, Snun/nbin = [6(3.23- l)]5a/ao, 
and the coefficient in brackets is 13.4. Both coefficients 
are large, implying the need for control of the fractional 
cross-section error at the 1% level for effective geometry 
determination. In the conventional approach that level 
of precision has not been achieved. 

Absence of extrapolation constraints can result in large 
systematic errors, as sketched in Fig. [TU The fractional 
error in n par t resulting from an incorrect endpoint value 
increases to 50% for peripheral collisions. In the conven- 
tional plotting format of the left panel (c/. Fig. 8 of (25j |) 
the impact of the systematic error (not uncertainty) is 
minimized by the nonlinear relation of n par t to fractional 
cross section. The 20% most peripheral collisions are 
omitted (below the left dashed line), and large system- 
atic error bars (uncertainty estimates) are applied be- 
low 50%. That strategy abandons all information on the 
transition from N-N to mid-peripheral A-A: how heavy 
ion collisions become distinct from elementary collisions. 
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FIG. 14: The fractional systematic error in n par t resulting 
from an incorrect lower endpoint (3/4 vs 1/2) plotted against 
ripart (left panel) and v (right panel) , with the same two values 
of the fractional cross section indicated in each panel. 

In contrast, the plotting format in the right panel pro- 
vides precise access to the peripheral transition region, 
provided systematic errors are brought under control. 

In the power-law context adequate control is accom- 
plished by invoking extrapolation constraints: aside from 
fluctuation effects both Glauber running integrals should 
extrapolate to endpoints ^1/2 to be consistent with the 
power-law distributions on n c h, pt and E t . On the other 
hand, both Glauber parameters should extrapolate to 1 
in the limit of N-N collisions or cr/co = 1. Parame- 
terizations satisfying those constraints reduce systematic 
uncertainties from the Glauber parameters to about 5% 
in the troublesome peripheral region (y ~ 2). 

As noted, systematic errors for path-length v are 
much reduced because of common-mode error reduc- 
tion. The relative error is 5v/v = {6 ■ 2.23/n^f — 4 • 
2.72/(n par t/2) 1/ ' 4 } 5a/ao. The coefficient in curly brack- 
ets has limiting values 2.5 (peripheral) and 1.25 (central). 
Thus, the error in v relative to the fractional cross sec- 
tion is 0(1-2) compared to relative errors for n part /2 and 
nun which are 0(10-15). With extrapolation constraints 
the systematic uncertainty in v can typically be reduced 
to < 2%, even in the peripheral region. 



C. Data errors 

Trigger and vertex-reconstruction inefficiencies, the 
latter typically significant for peripheral collisions, result 
in distortion of the minimum-bias distribution leading 
to systematic errors in the inferred centrality. The con- 
ventional method does not utilize critical a priori infor- 
mation available from p-p collisions which constrains the 
form of the minimum-bias distribution in the peripheral 
region. Without such constraints the fractional cross- 
section uncertainty can be 5-10% for peripheral collisions, 
making the n part and nun parameters meaningless in 
that region (50-100% error). For that reason the 20% 
most-peripheral part of the cross-section, which contains 
critical information on the transition from N-N to A-A 
collisions, is typically abandoned in the conventional ap- 



15 



proach. 

Trigger/ vertex systematic errors are illustrated in 
Fig. [T5] (left panel) . The overall event efficiency (trigger 
plus vertex reconstruction) is typically uniform and close 
to 100%, except for peripheral collisions where both effi- 
ciencies may be much reduced and strongly varying. The 
upper endpoint hq is precisely known (< 2%) from the 
power-law minimum-bias distribution (but not the con- 
ventional method). The lower endpoint n p ~ njvAr/2 can 
be inferred from independent experiments. Those two 
numbers provide a precise extrapolation reference for the 
power-law relation between n c h and cr/coi a $ illustrated 
in Fig. [15] (right panel) . 
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FIG. 15: Left panel: A cartoon of the power-law differential 
cross section illustrating typical trigger/vertex inefficiencies. 
Right panel: A cartoon of running integrals illustrating sys- 
tematic errors resulting from inefficiencies in the left panel. 



Assuming a modest centrality error Sa/ao ~ 1%, 
what is the corresponding relative error of the extrap- 
olated n^ 4 ? The relation is 5n p /n p = 45rip /rip — 
4(nJ /4 /np /4 - l)5a/a Q . For Au-Au at 200 GeV in two 
units of pseudorapidity Sn p /n p ~ 4((1050/4) 1/4 - 1) x 1% 
= 12%. By reverse argument we conclude that if one 
knows the N-N multiplicity to about 10% one can limit 
systematic error in the fractional cross section for data to 
f % for peripheral collisions, another extrapolation con- 
straint. That constraint in turn greatly reduces uncer- 
tainty in the Glauber parameters. The most peripheral 
bin is thereby restored to precision centrality studies. 



XI. DISCUSSION 

To unravel the complex dynamics of heavy ion col- 
lisions we must distinguish overlapping contributions to 
the hnal-state momentum distribution from initial parton 
scattering, subsequent in-medium parton dissipation and 
fragmentation and bulk-medium dynamics and fragmen- 
tation. Such a decomposition requires differential single- 
particle spectrum and two-particle correlation analysis 
precisely related to centrality over the entire centrality 
range from N-N to central A- A collisions. We require pre- 
cise knowledge of participant and binary-collision num- 
bers, nucleus overlap geometry and mean participant 



path length. In this study we have compared conven- 
tional and power-law centrality methods for Glauber 
model simulations, the Hijing Monte Carlo and RHIC 
data. 



A. Advantages of the power-law centrality method 

The power-law method offers the following advantages: 
f ) precise visual access to the entire minimum-bias distri- 
bution for any collision observable, 2) well-defined end- 
points n p and uq of the minimum-bias rich (or p t , E t ) 
distribution, with exact correspondence to the endpoints 
of the Glauber n part /2 distribution, leading to 3) a priori 
extrapolation constraints which improve centrality accu- 
racy for peripheral collisions up to f Ox , 4) simple and ac- 
curate paramcterizations of Glauber parameters, 5) flex- 
ible definition of centrality binning on n c h, and 6) useful 
extrapolation of the minimum-bias distribution to N-N 
collisions, even if the distribution is severely distorted by 
inefficiencies or backgrounds (next subsection). Those 
features provide first access to the most peripheral 20% 
of RHIC A-A collisions. 



B. Recovering from measurement distortions 

If the minimum-bias distribution is badly distorted by 
triggering/vertex inefficiencies or backgrounds the power- 
law method provides fall-back determination of collision 
centrality with good accuracy. We can reconstruct the 
minimum-bias distribution on n c h from a few pieces of 
information. Assume that f) the upper half of the distri- 
bution is undistorted except for smoothly-varying track- 
ing inefficiency, 2) the tracking inefficiency e for periph- 
eral collisions is known to f 0% and 3) the p-p multiplicity 
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is known to 5%. From the upper half of the distribu- 
tion we obtain rig and the mean slope of the distribution. 
From the p-p data and the measured efficiency we obtain 
n p = 0A5en pp and Uq = 0.45 en par tfi/2 ■ n pp . Those pa- 
rameters provide the linear participant-scaling reference. 
The systematic uncertainty in v given those conditions is 
< 2% for peripheral collisions. 

From n p and no we estimate the linear power-law ref- 
erence. Given the mean slope of the differential distribu- 
tion from its upper half and Eq. (|A3|) we can estimate the 
curved trajectory which lies between the linear limits. We 
can use either the resulting complete parameterization or 
the running integral of the data adjusted so that it agrees 
tangentially with the participant-scaling reference in the 
peripheral region to determine centrality based on n c h 
accurate to < 2% over the full centrality range. Exam- 
ples of such extrapolations and reconstructions are found 
in Figs. [31 [5] and [T3]of this paper. 
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C. Precision study of particle, p t and Et production 

Precise study of particle, p t and E t production per 
participant relative to centrality is essential to sepa- 
rate initial-state production mechanisms (parton scat- 
tering), in- medium parton dissipation, bulk- medium dy- 
namics and final-state hadron production (parton and 
bulk- medium fragmentation) . Collision dynamics should 
be confronted separately for total n c h, pt and E t in a large 
acceptance vs path length v. In contrast, (pt) or (et) 
(lower curves in Fig. [12]- right panel) measure p t and E t 
per final-state hadron, and are thus strongly modified by 
the final-state hadronization process, mixing initial-state 
and final-state collision mechanisms. 

We have shown in Figs. fTTT [T2l and IT31 that the power- 
law centrality method reveals new details of the colli- 
sion process at the few-percent level and reverses some 
apparent contradictions. We have related several mani- 
festations of the two-component model which provides a 
simple linear reference for production mechanisms. Sig- 
nificant deviations from that model, obtained for the first 
time with the power-law method, offer further details of 
the A-A collision process and bulk-medium properties. 
The per-participant distributions on v provide the means 
to look beyond hadronization into the dynamics of the 
pre-hadronic QCD medium. 



XII. SUMMARY 

We have introduced a new technique for determin- 
ing heavy ion collision centrality based on an observed 
power-law trend in the minimum-bias distribution on par- 

—3/4 

ticle multiplicity n c h of the form n h . The power-law 

trend implies that transformation to nj^ (power-law for- 
mat) should result in a nearly uniform distribution with 
physics-related variations confined to about ~ 20% of the 
mean value. The linear power-law plotting format pro- 
vides precise visual access to distribution structure. The 
clearly-identified endpoints (half-maximum points) of the 
minimum-bias data distribution correspond to the end- 
points of the participant-nucleon and binary-collision dis- 
tributions, providing a precise relation between data and 
Glauber simulations used to connect collision geometry 
to measured collision observables. 

We have confirmed that the power-law plotting format 
is applicable to Glauber nun and n part minimum-bias 
distributions as well as to data n c h, Pt and E t distribu- 
tions. We find that the minimum-bias distribution on 
n-part is almost exactly oc n pa 3 / t 4 , explaining the similar 
trend in data, and the distribution on nun is approx- 
Distributions on nj^ rt and n\j^ x are 
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imately oc n bin ' . 

therefore nearly uniform, providing simple and precise 
linear representations of the Glauber parameters vs frac- 
tional cross section. 

We have shown that application of power-law tech- 
niques to A-A centrality determination can reduce sys- 



tematic errors for peripheral collisions to the percent 
level, even when measurements are severely distorted 
by inefficiencies. The sharp reduction in systematic un- 
certainties results from application of extrapolation con- 
straints accessible only in the power-law context. 

We have applied the new centrality techniques and 
power-law parameterizations of the Glauber parameters 
to studies of particle, pt and E t production in Hijing 
and RHIC collisions. We have shown that the power- 
law method can reduce systematic errors in E t produc- 
tion of up to 50% to systematic uncertainties of about 
5%. We have also included three appendices which re- 
view particle-production algebra and numerical integra- 
tion techniques required to optimize the accuracy of col- 
lision geometry in the power-law context. 

TAT appreciates helpful discussions with J. C. Dun- 
lop (Brookhaven National Laboratories) and R. L. Ray 
(University of Texas, Austin) which led to Appendix IC 31 
This work was supported in part by the Office of Science 
of the U.S. DoE under grant DE-FG03-97ER41020. 



APPENDIX A: PARTICLE PRODUCTION 

According to a two-component model of particle pro- 
duction {n c h, as well as pt and E t production) centrality 
trends may follow a mixture of participant and binary- 
collision scaling [l7|. The two-component model is gov- 
erned by a single parameter x and expressed in the form 

2n c h/n part = n NN (l + x[u-l]), (Al) 

where njvjv is the N-N multiplicity and v is the mean 
participant path length [171 ] . The model has manifesta- 
tions in the differential cross section, its running integral 
and the per-participant-pair particle-production trend. 
We compare data to a participant-scaling reference in 
each case. We assume for simplicity that the participant 
minimum-bias distribution da/d(n part /2) 1 / 4 is uniform 
(no sinusoid component). Then, by the chain rule 
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where the first factor is the uniform participant-scaling 
reference, and the second factor is a Jacobian represent- 
ing the two-component model which we now calculate. 
Given the two-component relation Eq. (|A1[) we have 
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where for this derivation we have used the approximation 
v ~ (npart/2) 1 ' 3 (cf. Fig. [5]- left panel, dash-dot curve) 
consistent with the approximation nun ~ (npart /2) 4//3 . 
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The first two factors on the RHS form the participant- 
scaling reference. Eq. (|A3[) precisely describes centrality- 
determination plots for corrected data, such as Fig. [TOl 
(right panel). 

The required Jacobean factor is obtained from 

i/4 j. ~ /9W4 " (1 + -I--1]) 174 (A4) 
+ (r W /2) 1 / 4 (l + i[ v - l])- 3 / 4 W3. 

Introducing the Jacobian into Eq. (|A2j) we obtain 



1. Binning 

Confusion may arise in comparing the histogram of 
a density on a binned continuous variable with a dis- 
tribution on a discrete (integer) variable. A histogram 
bin entry represents the integral of a density over the 
bin. When divided by a bin width, the entry estimates a 
sample of the density within the bin. If a precise corre- 
spondence is sought between distributions on continuous 
and discrete variables (e.g., E t and n par t/2) or two dis- 
crete variables (e.g., n c h and n part /2), optimized binning 
and integration definitions are required. 



Kh n NN (Kar t , /2] 1/4 ~ KarUp/2] V<) 

(l + x[v- 1]) 3/4 

X 1 + x{v - 1) +2^5/4/3' 

which for x > predicts an approximately linear reduc- 
tion of the differential cross section on nj^ extending 

from n]^% to n^ 4 . Assuming x ~ 0.08 for 130 GeV 
RIIIC Au-Au data we predict a fractional reduction for 
central collisions of about 20% relative to the uniform 
participant-scaling reference (first factor on the RHS). In 
Fig. [3] we indeed observe a slope of ~ —0.15 with ~20% 
linear reduction over the full centrality range. 

We have thus related three examples of power-law cen- 
trality determination — the minimum-bias differential dis- 
tribution, its running integral and the per-participant- 
pair particle-production trend — with a system of model 
functions having two parameters: njvjv and x. Experi- 
mentally, the measured combination of tin n and no de- 
termines x. For instance, at 200 GeV no as 191 7ijviv(l + 
4.95 x). The x value so inferred must be consistent with 
the negative slope of the power-law minimum-bias dis- 
tribution, the curvature of its running integral and the 
positive slope of the particle-production trend. The last 
format, through the running integral, reduces the short- 
wavelength statistical noise on the minimum-bias dis- 
tribution at the expense of introducing long-wavelength 
noise from uncertainty in the Glauber n par t/2 parame- 
terization. The same relationships are true for p t and E t , 
as demonstrated in Figs. [T2l and [T3l ( although x is unique 
for each case). 



APPENDIX B: NUMERICAL INTEGRATION 

A major source of centrality error is relating Glauber 
parameters to measured quantities through running in- 
tegrals of the fractional cross section. We want to reduce 
those errors by improving the running-integration meth- 
ods. To that end we review some details of binning and 
numerical integration. 



2. Power-law transformation 

The elements of a minimum-bias distribution are event 
counts Ni, where i is a bin index on continuous variables 
p t or Et or labels values of discrete variables n c h, n part /2 
or ni,in- For bin widths Sxi on x the minimum-bias 
density is estimated by dNi/dx = Ni/Sxi. The power- 
law form is obtained by the transformation dN/dx 1 ^ 4 = 
4x 3 / 4 dN/dx. The power-law distribution on multiplicity 
is dN/dnJ 4 = AnJ 4 dN/dn ch = An J 4 N rich . The conven- 
tional minimum-bias distribution is a non-uniform distri- 
bution on a uniform bin system. The power-law form is 
a nearly-uniform distribution on a non-uniform bin sys- 
tem. The latter makes details of the density distribution 
and its integration more accessible and facilitates extrap- 
olation. 



3. Rebinning 

When transforming a histogram from space x to space 
y rebinning may be desirable. Assume a distribution fi 
(e.g., event number) on discrete variable n.j. We want a 
binned distribution on continuous variable y with uni- 
form bin widths 8y, transformation yi — n 1 ^ 4 and values 

<7, = An'J 4 fi . We define uniformly-spaced bin edges jjk 
on y with index k. For rebinning we step through values 
ni and histogram elements gi. Starting with bin k = 1 on 
y we sum histogram elements gi within the k th bin into 
Gfc, number of steps on n into and values rii into N/~ 
(this could also be values n~J 4 ). We test the y bin edge 

for each nf. if n 1 ^ 4 < we continue the n loop. If not, 
we advance to y bin k + 1 and continue the n loop. We 
increment k until the end of the specified y interval. We 
then form Gk/Mk = y~k and N^/Mk = fik (this could also 
be n 1 / 4 ^) as uniformly-spaced transformed densities gk 

and means n]/ 4 (or n 1 / 4 ^) on bin centers yk = yk~ Sy/2. 
The minimum-bias distribution {fi} on particle multi- 
plicity rich is thus rebinned to a 'power-law' distribution 
{gk} on n 4 J 4 , as in Fig. O (left panel). 
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4. Numerical integration 

To form the Riemann sum of function f{x) over a 
closed interval [x a ,i£b] of continuous variable x the in- 
terval is partitioned into M bins i of width 5xi (not nec- 
essarily equal), and sample points chosen within 
the bins to sample the density f(x) — ► fi = f(xi). The 
Riemann sum is F(x a ,Xb) = F m = J2iLi^ x i hi where 
x a is the lower edge of the first bin and Xb is the upper 
edge of the last bin. Depending on the choice of sample 
positions within the bins the Riemann sum definition is 
continuously variable from 'lower' to 'upper' sum. Typ- 
ically, the X% 3X6 bin centers ('middle' Riemann sum), 
but we are free to adjust sample positions relative to bin 
edges to insure mutual compatibility of running integrals. 

The Riemann sum on discrete variable n of distribu- 
tion fi at integer values rij with i £ is apparently 
simple: Fm = J2i=i fi- However, if we extend n to a con- 
tinuous variable with unit-width bins centered on integer 
values the Riemann sum can be adjusted by shifting the 
bin edges relative to fixed integer positions (as opposed 
to shifting sample points within fixed bins). We require 
that flexibility to insure compatibility of various running 
integrals defined below. 



5. Running integrals 

To apply the Glauber model precisely we must define 
the running integral on integer variables (e.g., n part /2) 
consistently with integrals on other integer variables 
(e.g., n c h) or continuous variables (e.g., pt or E t ). To 
insure compatibility among running integrals on discrete 
and continuous variables we treat discrete variables as 
continuous binned variables. The running integral on a 
continuous variable is F(x a ,x m ) — F m = X)i=i $ x i /*> 
where x m is the upper edge of the m th bin and m £ 
[1, M]. The running integral of the distribution must be 
properly related to its independent variable. The corre- 
sponding running integral of the independent variable is 
x m = x a + J27Li & x ii assuming that integration is be- 
tween leading and trailing bin edges. 

If / on n is transformed to f 1 ! 4 on n 1 ^ 4 the bin 
edges are not symmetrically placed about the sample 
point, and individual bin edges and bin widths must 
be determined explicitly to minimize discretization er- 
rors. We use parameters u and v with u + v = 1 to de- 
note bin-edge positions relative to integer positions on 
n. Thus, bin edges relative to m are at n,; + u and 
rii — v. From the bin edges we calculate bin widths 
Sn 1 / 4 = (n + it) 1 / 4 — (n — v) 1 / 4 . (u,v) can be ad- 
justed to satisfy extrapolation constraints. E.g., u = 0.45 
rather than 0.5 due to the NBD distribution skewness and 



= + u = 0.45. The running integrals are then 

m 

nil 4 = nl/ 4 + J2Snl /4 = (n m + U y/ 4 (Bl) 
i=l 

m 

F m = £// /4 -{K+^) 1/4 -K- V ) 1/4 }- 

»=i 

We plot the running integral and running limit as F m vs 
nU, 4 (e.g., Fig. [6|. It should be clearly indicated whether 
a plot represents bin entries and bin centers or running- 
integral sums and corresponding bin edges. Because of 
the strong nonlinearity of n 1 / 4 , exact bin-edge positions 
relative to the integers are very significant for small n 
and irrelevant for large n. 

APPENDIX C: RELATING CENTRALITY 
PARAMETERS 

The underlying premise of centrality determination is 
that collision observables and Glauber geometry parame- 
ters are related statistically by joint probability distribu- 
tions which are not directly observable. What are acces- 
sible are projections of the joint distributions onto their 
margins. Pairs of variables (measured and Glauber) can 
then be related by running integrals of the marginal dis- 
tributions. We now discuss technical details of that pro- 
cedure. 



1. Joint distributions 

The basis for relations between n c h, Pt and E t on the 
one hand and n part /2 and nu n on the other is joint den- 
sity distributions on pairs of variables: a measured quan- 
tity and a collision-geometry parameter (e.g., n c h and 
n part / 2 ). The locus of means (curve describing condi- 
tional means) provides the functional relation between 
two quantities. Joint distributions are not directly ob- 
servable experimentally, but one projection (marginal 
distribution on n c h, Pt or E t ) is measured, and the other 
(marginal on n part /2 or Ubin) is obtained from a Glauber 
Monte Carlo simulation. We relate one variable to the 
other through running integrals of the marginal projec- 
tions. Each normalized running integral is the fractional 
cross section in the form 1 — ct/ctq- In the power-law con- 
text precise access to the locus of means is provided by 
the endpoints or half-maximum points of the marginal 
(minimum-bias) distributions. 

For joint distribution p(x, n) there are two conditional 
distributions: (x) n (mean of x for fixed n) and (n(x)} 
(mean of n for fixed x). We expect (x) n and (n(x)) 
to correspond for intermediate values of x and n, but 
to diverge near the joint distribution endpoints due to 
'fluctuations' (finite width of the joint distribution). In 
the power-law context the nearly-linear relation between 
centrality parameters over most of the joint distribution 
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can be determined precisely by the marginal endpoints 
determined as the half-maximum points of the marginal 
distributions. 



2. n part scaling and n ch 

A key element of A-A centrality determination is the 
participant-scaling (wounded- nucleon) hypothesis: a pro- 
jectile nucleon can 'participate' in at most one N-N colli- 
sion, and each such participant pair produces final-state 
hadrons as in an isolated N-N collision. In that limit 
the A-A multiplicity 'scales' as n c h oc n par t/2. Alterna- 
tively, in the binary-collision limit A-A multiplicity scales 
as rich oc nun- For strict participant scaling n part can- 
not have odd values; e.g., there cannot be exactly three 
participants in an A-A collision. Once a nucleon has par- 
ticipated in an N-N interaction (is 'wounded') it cannot 
produce a third participant in a second encounter, even 
though the number of binary collisions does increase by 
one in the binary-collision context. In the participant- 
scaling limit rip ar t/2 should assume only integer values, 
and in this paper we treat that symbol combination as 
the basic statistical variable. 

We can combine participant scaling with the power-law 
form of the minimum-bias distribution to define a precise 
relation between n c h and n part /2 valid to < 5% for pe- 
ripheral collisions. In the limiting case n c h oc n par t/2 
is strictly obeyed (ignoring N-N fluctuations). We em- 
ploy an extrapolation constraint at the lower endpoints 
to provide precise registration of n c h and n par t/2, even 
if the measured minimum-bias distribution on n c h is 
severely distorted by inefficiencies and backgrounds. If 
N-N fluctuations are symmetric about the mean then the 
minimum-bias distributions on n par t/2 and n c h extrap- 
olate to n parttP /2 — 1/2 and n p = 1/2 • n^pf respec- 
tively. The upper endpoints of the participant-scaling 
reference go to n par tfi/2 and no = n part .a/2 ■ njvAr, where 
n pa rt,o/2 ~ 191 at yfs — 200 GeV. In the next section we 
justify those endpoint values with simulations and show 
how to achieve precise correspondence for real N-N fluc- 
tuations with nonzero skewness. 



3. Running integrals and N-N fluctuations 

N-N fluctuations can change the position of endpoint 
n p relative to tinn, or equivalently the value of n par t tP . 
To understand the systematics we model the joint dis- 
tribution on (n part ,n c h). We assume strict participant 
scaling. We also assume that the differential cross sec- 
tion is constant on n part /2 rather than (npart/2) 1 / 4 . The 
projection onto n par t/2 is a set of delta functions on the 
integers. The projection onto rich depends on N-N fluc- 
tuations, as shown in Figs. [16] and [17] We relate n c h to 
npart/2 through running integrals of the projections. 

In Fig. [16] (upper-left panel) we show the projected 
density on n c h for a narrow gaussian N-N multiplicity dis- 
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FIG. 16: Upper-left panel: Minimum-bias distribution on n c h , 
with narrow gaussian N-N distribution. Upper-right panel: 
running integral of the minimum-bias distribution on n c h to 
obtain the fractional cross section. Lower-left panel: Con- 
ditional distribution of (n par t/2) on fractional cross section. 
Lower-right panel: Conditional distribution of (n part /2) on 
n c h- The open circles denote the complementary conditional 
distribution (n c h) on n part /2. 



tribution and 1-5 participant pairs, assuming umn = 5 
and a = 1 < ^tinn- In the upper-right panel we show 
the running integral of fractional cross section 1 — a/ao on 
rich- In the lower-left panel we show the conditional dis- 
tribution (n part /2) on a /ao (different from the running 
integral on n par t/2), and in the lower-right panel we show 
the conditional distribution (n par t/2) on n c h- The open 
circles indicate the complementary discrete conditional 
distribution (n c h) on n part /2. The two are consistent 
by construction. The lower-right panel reflects a 'chain- 
rule' relating (n par t/2) to n c h through the lower- left and 
upper-right panels. 

In Fig. [T7] we show the same set of plots for a = 3 > 
y/ n NN- The smooth variations allow us to examine the 
structure at the endpoints. In the upper-left panel we 
see that the lower half-maximum (dotted line) falls at 
n p = njvw/2 = 2.5 (n part p — 1/2). That endpoint is 
stable over a range of N-N widths for symmetric N-N 
fluctuation distributions. In the upper-right panel the 
running integral endpoints (extrapolated values) corre- 
spond to the half- maximum points of the differential dis- 
tribution. In the lower-left panel the lower endpoint is 
(n P art,p/2) = 1/2 as expected. 

In the lower-right panel the locus of means for the two 
complementary conditional distributions (solid curve and 
dots) agrees over the central region, as determined by the 
match of endpoint values in the lower-left and upper-right 
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FIG. 17: The same plots as in Fig. \W\ but with broader N- 
N multiplicity distribution to illustrate asymptotic behavior 
near the endpoints for RHIC data. 

panels. For real systems we must adjust the Riemann 
sum definitions so that (for instance) the locus-of-means 
endpoints on n c /j and n part /2 coincide. Also, we must 
define and integrate n part /2 and ribi n so that both lower 
limits are the same to insure that v goes to 1 in the limit. 

The minimum-bias distribution for real data is based 
on an N-N multiplicity distribution with nonzero skew- 
ness, commonly modeled by the negative binomial distri- 
bution (NBD). The NBD is defined by 

P(n; ^ )= |^il (p - 1)y ' (C1) 

with fixed k. For comparison, the binomial distribution 
with the same notation is 

p{n ' N ' p) = w^ {p ~ l)kpn (C2) 

with fixed TV = k + n. The NBD has mean fj, = kp/(l—p) 
and variance a 2 = /^/(l — p), whereas the binomial has 
H = Np and variance a 2 = p(l—p). Both go to a Poisson 
distribution with o 2 = p if N, k — > oo and p — > with 
fixed n. For the NBD (a 2 NBD - a 2 Pmsson )/n = \ijk = 
p/(l—p) is a measure of normalized variance (fluctuation) 
excess relative to Poisson due to correlations. For N- 
N collisions in the STAR TPC acceptance /z ~ 5 and 
fc~ 10, withp~ 1/3. 



We want to understand the effect of skewness and 
width changes on the value of n partiP . The NBD skew- 
ness depends on its control parameter k. If pb/k — *■ 
the NBD goes to a Poisson distribution. For larger p/k 
the NBD is increasingly skewed and its width increases. 
The lower half-maximum point on the differential cross 
section then moves below tinn/2 and the participant- 
pair-number running integral definition must change to 
accommodate. 




FIG. 18: Model of a minimum-bias distribution on n c h with 
NBD N-N distribution, showing the changing position of end- 
point n p ~ 2.5 relative to tijvjv = 5 as the NBD correlation 
parameter k is varied. For RHIC data (k ~ 10) n p ~ 0.45 umn 
is a good approximation. 

In Fig. [18] we compare minimum-bias distributions 
based on NBD and gaussian distributions. In the left 
panel we show a minimum-bias distribution (solid curve) 
with the individual N-N distributions modeled by an 
NBD with k = 100 and unn — 5 (dashed curve) . For 
comparison the dotted and dash-dot curves are based on 
a gaussian N-N model, as in the previous figure but with 
a = V5 = 2.25. The value of n p for the NBD is slightly 
larger than the njvjv/2 = 2.5 for the gaussian. The NBD 
with k = 100 is nearly Poisson, with a only slightly larger 
than y/Jl. 

In Fig. [18] (right panel) we show similar distributions 
but with k = 10 corresponding to data. Due to the in- 
creased skewness and width of the NBD with the smaller 
k the value of n p is shifted below tlnn/2 to n p « 
0.45 tinn ■ We conclude that n part , p /2 — 0.45 ± 0.05 is a 
good match to RHIC A-A data. We therefore shift bin 
edges on n part /2 down slightly relative to the middle Rie- 
mann sum. Equivalcntly, we adjust the n part p /2 value in 
the parameterization of Eq. ^ and adopt n p « 0.45njv7v 
as the extrapolation constraint for minimum-bias distri- 
butions on n c h- The 10% relative uncertainty in n par t, p /2 
propagates to a mere 0.7% uncertainty in a/ao for pe- 
ripheral collisions. However, it corresponds in particle, 
Pt or E t production plots such as Fig. [13] to a 5% uncer- 
tainty localized near v = 2. 
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